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Studying the leptonic structure of galaxy cluster atmospheres 
from the spectral properties of the SZ effect 
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Abstract. We study the energetics of galaxy cluster atmospheres by analyzing the SZ effect spectra around the 
crossover frequency. We calculated analytically the expressions of both the crossover frequency and the spectral 
slope of the SZE around the crossover frequency in various cases: a thermal electron population; a power-law, 
non-thermal, electron population; and a population of electrons experiencing a stochastic acceleration. We find 
that the value of the crossover frequency Xo of the SZE depends significantly on the cluster peculiar velocity V r 
which determines the amplitude of the kinematic SZE), while the value of the slope of the SZE does not depend 
on the kinematic SZE spectrum in the optimal frequency range around the crossover frequency of the thermal 
SZE, i.e. in the frequency range x = 3.5 — 4.5. Therefore, while the amplitude of the SZkin produces a systematic 
bias in the position of the crossover frequency Xo, it does not affect significantly the spectral slope of the SZE. 
We therefore propose using measurements of the spectral slope of the SZE to obtain unbiased information about 
the specific properties of various electron distributions in galaxy clusters as well as in other cosmic structures in 
which a SZE can be produced. 
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1. Introduction 

The SZ effect (hereafter SZE) is a promising tool for 
studying the complex physics of cluster atmospheres 
because it is sensitive to the specific features of the 
spectra of the various electron populations responsi- 
ble for the CMB photon Comptonization (see e.g. 



plasmas. 

For the majority of galaxy clusters the most prominent 
form of SZE is the thermal one, which provides a CMB 
intensity change of 



th 



Colafrancesco (2007) for a review). A SZE can be pro- 
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duced, in fact, by different types of leptonic plasmas 
in astrophysical environments, such as: thermal (hot 
and/or warm) electrons in the atmospheres of galaxies 
and galaxy clusters ( Birkinshaw (1999)] Itoh et al. (1998)| 



Colafrancesco et al. (2003)); non-thermal (and relativis- 



as produced by Inverse Compton Scattering (ICS) of 
CMB photons off the thermal electron population re- 
siding in the cluster atmosphere with a Comptonization 
parameter y t h = (cr/meC 2 ) J d£n e kBT e given in terms 
of the electron plasma number density n e and temper- 



tic) electrons in clusters and in the cavities pro- 
duced by AGN radio lo b es pnsslin fc Kaiser (2000) [ 



Colafrancesco et al. (2003) Colafrancesco (2005) ); sec 



ondary electrons produced by Dark Matter annihilation 



in cosmic structures ( Colafrancesco (2004) I , in addition 



to the kinematic SZE related to the bulk motion of these 
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ature T e (see e.g. Birkinshaw (1999) for a review and 
Colafrancesco et al. (2003) for a general derivation) and 
in terms of fundamental constants, namely the Thomson 
cross-section ar, the electron mass m e , the speed of light 
c, the Boltzmann constant krj and the Planck constant h. 
The function g(x) depends on the adimcnsional frequency 
x = hv/k B T , where T = 2.726 K is the CMB temper- 
ature, containing all the spectral information about the 
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SZE (see e.g. Colafrancesco et al. (2003) for details) 



non-thermal SZE due to a relativistic plasma, the DM- 



Thcre are three basic spectral features that characterize 
the thermal SZE signal: 

i) a minimum in its intensity located at a frequency 

x th ,min « 2.265(1 - O.O9270 e + 2.380*) 

+t(-O.OO674 + O.4660 6 ) , (2) 

where e = kT e /m e c 2 , and r = ot J d£n e is the optical 
depth of the electron plasma, whose exact position de- 
pends weakly on the electron spectrum (i.e. on the elec- 
tron temperature T e and number density n e ) and equals 
- 2.26; 

ii) a crossover frequency, Xq, whose value depends on the 
electron pressure/energy density and the electron optical 
depth 



X th , « a(T e ) + Tb(T e ) 



(3) 



with a{T e ) = 3.830(1 + 1.1620 e - 0.81440;?) and b(T e ) = 

3.O210 e — 8.6720^, and is found at a frequency > 3.83 for 

increasing values of T e (the value X t h,o = 3.83 is found in 

the non-relativistic limit, or in the limit T e — > 0); 

iii) a maximum of its intensity whose frequency location 

is 



xth.max ~ 6.511(1 + 2.410 e -4.960 e 2 ) 

+r(O.O161 + 8.160 e - 35.90^) 



(4) 



and depends sensitively on the nature of the electron 
population and on its energy (momentum) spectrum (see 



Dolgov et al. (2001 for the case of electrons with a ther- 



mal spectrum; see also Colafrancesco et al. (2003) 


and 


Colafrancesco (2004) 


Colafrancesco (2005)| 


Colafrancesco (2007) for the case 


of electron popula- 



tions with different spectra). 

Since the frequency location of the zero in every SZE 
signal contains a crucial dependence on the pressure (or 
energy density) of the electron population (see e.g. Eq. [3] 
for the thermal case; see also ? for the discussion of the 
general cases), it has been widely proposed to use this 
property to measure e.g. the cluster temperature (for a 
thermal electron population) or the pressure/energy den- 
sity in the case of a more general electron population. 

However, observations of the SZE close to the crossover 
frequency are biased due to the presence of the kinematic 
SZE 



Ah 
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(see Sunyaev & Zel'dovich (1980) see also Eq. [8] below) 



which is associated with the peculiar velocity V r of the 
cluster. The unknown value of V r therefore limits the abil- 
ity to measure the cluster temperature (or more gener- 
ally its energy density) directly through the displacement 
of the crossover frequency Xq of the SZE in the correct 
relativistic treatment. The possible presence of multiple 
sources of SZE in galaxy clusters (e.g. the thermal SZE, 
the SZE due to an additional warm gas component, the 



induced SZE; see Colafrancesco (2007) for a recent review) 
could provide further uncertainties and biases in the mea- 
surement of the crossover frequency Xq of the dominant 
thermal SZE, since they add to the thermal SZE with spe- 
cific spectral shapes and amplitudes of their signals. 
This parameter degeneracy might be broken - in princi- 
ple - by observing the SZE signal at many frequencies, a 
procedure that is limited, however, by the observational 
setup of the current SZ experiments operating in a limited 
number (usually ~ 3-4) of frequency bands of relatively 
wide frequency width (~ 15 — 30% of the band central 
frequency) requiring precise inter-calibration. 

In this paper we explore the possibility of extract- 
ing unbiased information about the nature of the electron 
population that produces the SZE from the spectral shape 
of the relative SZE spectrum. 

One of the frequency region where large variation of 
the SZE spectra are found is around the crossover fre- 
quency of the thermal SZE, at frequency ~ 220 GHz (see 



Fig.l in Colafrancesco (2007)). Therefore, this is a favor- 
able region for performing a spectral analysis of various 
SZE components with the aim to separate the contribu- 
tions of various electron populations. 
To provide quantitative estimates of our proposal that 
could be tested by SZE experiments with spectroscopic 
capabilities, we will present an analytical derivation of 
the spectral slope of various sources of SZE around the 
crossover frequency of the thermal SZE, i.e. in the range 
x = 3.5 — 4.5 (corresponding to the frequency range 
~ 200 — 255 GHz). We show that, in this frequency re- 
gion, the SZE spectral slopes are mostly independent of 
the cluster peculiar velocity, at variance with the value 
of the crossover frequency Xq, that is inevitably biased 
by the unknown value of the cluster peculiar velocity V r . 
For the sake of illustration, we first discuss in Sect. 2 the 
cases of the non-relativistic and relativistic formulations 
of the thermal SZ effect. In Sect. 3 we discuss the case 
of a completely non-thermal, relativistic SZE for a power- 
law electron spectrum. In Sect. 4 we discuss the physical 
case of the SZE produced by a Maxwellian electron distri- 
bution experiencing a stochastic acceleration process. We 
discuss our results and draw our conclusions in Sect. 5. 
Throughout the paper, we use a flat, vacuum-dominated, 
CDM cosmology with H o =70 km s _1 Mpc -1 , fi M =0.3, 
Ov=0.7. 

2. The spectrum of the thermal SZE around the 
crossover frequency 

The spectra of the thermal SZE in the relativistic and 
non-relativistic limits differ quite significantly around 
the crossover frequency for increasing values of thermal 



plasma temperature (see e.g. Birkinshaw (1999) and ref- 
erences therein). This fact has been widely indicated as 
a possible means of extracting physical information, such 
as e.g. the electronic plasma temperature T e from a mea- 
surement of the crossover frequency X (T e ) of the SZE 
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spectrum, if the cluster peculiar velocity V r is known to a 
sufficiently high accuracy. However, cluster peculiar veloc- 
ity information is often uncertain or even missing, which 
means that this technique can be affected by strong sys- 
tematic biases. 

In the following, we show that measuring the spectral 
slope of the SZE spectrum around the crossover frequency 
provides a means of extracting information on the electron 
distribution that is independent of the cluster peculiar ve- 
locity. 

We compare the spectral slopes of the thermal SZE, 
SZ t h , and the kinematic SZE, SZ^m , calculated for an 
isothermal cluster in the frequency range 200-255 GHz, 
corresponding to the range x ~ 3.5 — 4.5 of the adimen- 
sional frequency x = hv/kT$. 

We begin, for the sake of clarity, our discussion from 
the non-relativistic case where the Kompaneets (1957) 
equation can be solved analytically and leads, following 
Eq. ([1]), to an analytic expression for the CMB intensity 
change in the non-relativistic regime 



with 

9nr{x) 



2(fcr ) 

(he) 2 



(e x - If 



kT P 



m e c 



'9nr (x) 



1 



-4 



(0) 



(7) 



The subscript 'nr' denotes the fact that the previous ex- 
pression was obtained in the non-relativistic limit. 

The kinematic SZE intensity change (Sunyaev & 
Zeldovich 1980), due to the bulk motion of the cluster 
gas with peculiar velocity V r along the los, writes as 



AI k = 
with 

hnr (#) 



(e x - l) 2 



(8) 



(9) 



In the following, we use the quantity i(x) to represent the 
adimensional CMB intensity spectral distortion 



hh,k{x) 



AI, 



th,k 



2(fcT ) 
(he) 2 



(10) 



which contains relevant spectral information on the ther- 
mal and kinematic SZE. 

The total SZE intensity in the case of a cluster with a 
single thermal electron population moving along the line of 
sight with a peculiar velocity V r is due to the combination 
of SZth and SZfcj„ and is given by 



i(x) 



V r m e c 2 ' 
SW U7^ h ( x ) 



(11) 



The series expansions of g(x) and h(x) around the adi- 
mensional frequency x — 4 (see Appendix A) provides the 




Fig. 1. The functions g nr (x) (solid) and h nr (x) (dashed) 
are shown in the frequency range x = 3.5 — 4.5. 



values of the spectral slope S and of the crossover fre- 
quency X . 

The slope of the SZE spectrum around the crossover fre- 
quency Xo is defined by the quantity 

i(x) - i(X ) _ i(x) 



S = 



(12) 



X — Xq X — Xq 

In Fig. 1, we compare the functions g(x) and h(x) in 
the frequency range of our interest. Since the SZkin spec- 
tral shape h(x) is almost fiat in the frequency range 
x ~ 3.5 — 4.5, the slope S of the thermal SZE can be 
used to extract unbiased information about the nature of 
the electron population, i.e., in this case, its temperature 
T e . 

In the non-relativistic limit, we find (see the Appendix 
for a derivation) that the slope of the total SZE (thermal 
plus kinematic) is given by 



kT e g(X + 5x) 



Ul e C 



5x 



4.25 % = 0.08 



m e c^ 



lOkeV 



(13) 



and does not depend on the cluster peculiar velocity but 
only on the electron temperature T e . 
The crossover frequency of the total SZE 

V r 10keV\ (w) 



X, 



O.total.nr 



= 3.83 + 0.193 



10 3 fcr 



obtained from the condition i(X ) = (see Appendix A.l) 
depends, instead, on both kT e and V r , so that a measure of 
the crossover frequency for the thermal SZE cannot pro- 
vide an unambiguous estimate of the cluster temperature 
T e . Specifically, for values kT e = 10 kcV and V r = ±1000 
km s _1 , the values of Xo ito tai,nr is found to be in the 
range 3.64 — 4.02, while for V r = one recovers the non- 
relativistic value X 0i totai,nr = 3.83, which, obviously, does 
not depends on the cluster temperature. 

In the relativistic treatment, the CMB spectral distor- 
tion due to the thermal SZE is given by the expression 



AL 



Hi 



(he) 2 



dsPi(s) 



cxp(xe 



1(15) 
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where 



P 1 (s)= / d0f e (j3)P(a,0) 



(16) 



with / e (/3) being the velocity spectrum of the electron pop- 
ulation and P(s, (3)ds being the probability that a single 
scattering of a CMB photon off an electron with speed f3c 
causes a logarithmic frequency shift s = ln{v' ' jv). Here, 
wc use the first order approximation in r of the pho- 
ton redistribution function P(s) for the sake of simplic- 
ity. Analogous conclusions hold, nonetheless, for the more 
general derivation (see e.g. Colafrancesco et al. (2003)). 



The series expansion of i{x) around the crossover fre- 
quency X was calculated in Eq. (A. 11) of the Appendix. 
Using the expression for the crossover frequency Xo(T e ) = 
3.830(1 + 1. 1626> e -0. 81446^) given by polgov et al. (2001[ 



one can calculate the slope of the thermal SZE in the rel- 
ativistic treatment as a function of kT e . A good approxi- 
mation of the slope S in the relativistic treatment, and in 
the temperature range kT e = 1 — 20 keV, is given by 



4.25- 



= 0.08- 



kT e 
m e c 2 
kT e 
WkeV 



-30 



kT e 
m e c 2 



1 - 0.138- 



kT e 
WkeV 



(17) 



and, again, depends only on the cluster temperature T e . 
The slopes of the non-relativistic and relativistically cor- 
rect thermal SZE spectrum are shown in Fig. [5] for com- 
parison. The value of the crossover frequency of the total 



0.15 



H 




0.05 



Fig. 2. The spectral slopes of the non-relativistic SZE 
(solid curve) and of the relativistically correct SZE spec- 
trum (dashed curve) calculated in the range x = 3.5 — 4.5 
for a thermal electron population are shown as a function 
of the electronic temperature T e given in keV. 

SZE in the relativistic treatment is given by 
V r h[X (T e )} 



Xo,total,rel — Xo(T e ) + 

« A (T e ) 



Srel 

0.193 



1-0.138- 



WkeV 




-0.5 



Fig. 3. The spectral shape of the normalized CMB inten- 
sity change i(x) evaluated for a power-law electron spec- 
trum with a minimum momentum p\ = 1 (dotted line), 3 
(dashed) and 10 (solid). 



V r 



I0 3 kms 



WkeV 
kT„ 



(18) 



(see the Appendix for a derivation) where Xo(T e ) = 
3.83(1+1.162-^-0.814-^). These results show that, while 
the value of the crossover frequency XQ yto tai,rei depends on 
both kT e and on V r (see Eqs. IT8l and fA.13j) . the value of 
S does not depend on V r in the frequency range around 
the crossover frequency Xq (see Appendix A for details). 
This is because the slope of the kinetic SZE is approxi- 
mately zero in this frequency range. Therefore, a measure 
of X for the thermal SZE cannot provide an unambiguous 
estimate of the cluster temperature. In contrast, a mea- 
sure of the slope of the thermal SZE can provide unbiased 
constraints on the value of the electron temperature and, 
hence, the nature of the electron plasma. 
In the following sections, we demonstrate that this is a 
general result for the SZE that can be applied to the case 
of other electron plasmas residing in cluster atmospheres. 

3. The spectrum of the non-thermal SZE 

We consider a non-thermal electron population with a 
power-law spectrum of the form 



fe( P ) = 
feip) = 



l-o 

Pi 







Pi <P<P2 

elsewhere 



(19) 



where f e (p) is the momentum distribution, and the mo- 
mentum p = j3j is normalized to unity. This electron spec- 
trum is expressed as a function of /3 by 



1 



1 



l-(3 2 



3/2 



.(20) 



The spectral shape of the quantity i(x) evaluated for the 
previous power-law electron spectrum is shown in Fig. [3] 
for different values of the minimum momentum p\ of the 
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Fig. 4. The dependence of -Xo,non-th on the lower cutoff 
momentum p\. We do not consider here the effect of a 
possible kinematic SZ effect. 

electron spectrum with a value of the momentum p2 = 
1000. 

The approximate value of the crossover frequency of the 
non-thermal SZE for a monoenergetic electron spectrum 
in the limit 7> 1 has been derived in the Appendix, and 
is given by 



L 0,non— th 



2ln 



(21) 



The dependence of Xo tnon - t h on the lower momentum p\ 
of the power-law electron spectrum is shown in Fig. [4] (this 
figure show the values of Xo,non-f/i without the effect of 
the peculiar velocity). 

The quantity Xo jnm -tii is assumed to have values 
larger than x ?J 5 for all values p\ > 1 (see Fig. @|. This 
means that for a system containing only a non-thermal 
electron population, the value of Xo. non _ t h (and hence of 
the energetics/pressure of the relative electron population) 
cannot be determined without a precise knowledge of its 
peculiar velocity V r . 

The value of the slope of the non-thermal SZE is also 
affected by the presence of a kinematic SZE in this case. 
However, both the SZj. in and primary CMB anisotropics 
have a flat spectral shape in the frequency range x = 
3.5 — 4.5, so that their contribution to the slope of the 
total SZE is negligible in this frequency range. 

3.1. Combination of thermal and non-thermal SZE 

For the case of realistic clusters with radio halos and/or 
X-ray cavities (i.e. those clusters definitely containing a 
population of relativistic electrons whose energy spec- 
trum is often represented by a power-law distribution, 
see Colafrancesco (2007) for a review), there is usually a 



cospatial distribution of thermal and non-thermal plas- 
mas. We calculate the values of the cross-over frequency 
Xq and the slope 5 of the SZE for a combination of a ther- 
mal plus a power-law electron spectra with the addition 
of a peculiar velocity of the cluster. 



In this case, we found that the value of the crossover 
frequency Xq depends on the cluster temperature, the pe- 
culiar velocity, and the optical depths of the thermal and 
non-thermal electron populations 



X = 3.83 + 0.193 



Vr 



WkeV 



+0.148 



10 3 km/s kT e 
T re i/r 10keV^ 

0.01 W e 



(22) 



(sec Appendix for details). 

The value of the slope of the total (thermal plus non- 
thermal plus kinetic) SZE is given by 



S = 0.084. 



1 + 7-10 



-3 



Vr 



10keV 



+0.038 



lOkeV 
T re i/r lOkeV 



10 3 km/s kT e 



0.01 kT e 



(23) 



To estimate specific values of Xq and S for a repre- 
sentative Coma-like cluster that contains both thermal 
and non-thermal plasmas, we choose representative val- 
ues kT e = 8.2 kcV, pi = 10, and T re i/r = 0.01. 
For the single Maxwellian spectrum, one obtains in this 
case the reference values Xq = 3.9 and S = 0.058. The 
zero of the total SZE evaluated for a combination of ther- 
mal and non-thermal populations is Xq = 4.09, and the 
increase in the value of X is in agreement with Eq. (|22l) . 
The spectral slope of the total SZE for a combination of 
the two electron populations is S = 0.061, and the mild 
increase in S in this case is also in good agreement with 
Eq. CU). 

The contribution of the kinematic SZE to the deter- 
mination of Xq cannot be ignored. Using Eq. (|2"2")) . we cal- 
culated the contributions to the value of X of the terms 
depending on V r and r rei / r . The contribution of V r to the 
value of Xq can be neglected in contrast to the contribu- 
tion of the non-thermal electron population only for val- 
ues \V r \ < 776.8kms- 1 (^f^); therefore, even for values 
(Trei/r) ~ 0.01, the value of Xq is substantially affected 
by possible values of V r in the range ~ 800 — 1000 km/s, 
as found in the tail of the peculiar velocity distribution 
(with rms value of « 300 ± 80 km/s) of galaxy clusters 



(see Giovanelli et al. (1998) 



In contrast, (from Eq. I23[) . the impact of the kinematic 
SZE (i.e. of the cluster peculiar velocity) on the slope 
can be neglected in comparison with the impact of the 
non-thermal electron population, for values | | < 5.4 • 
10 3 kms~ 1 ( T q e q( t ) ; this means that the slope of the total 
SZE is almost unaffected by realistic values of the cluster 
peculiar velocity which are well below 5000 km/s. 

4. A specific test case: the SZE spectrum of 
accelerated electrons 

To test the power of the method discussed in this paper, 
we apply our previous analysis to the specific case of a 
cluster in which the SZE signal is produced by an electron 
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population experiencing a stochastic acceleration process 



from a background plasma (see Dogiel et al. (2007) for an 



extensive discussion of this scenario). 

To describe the distribution function of the acceler- 
ated electrons, we use a kinetic equation modeling the 
influence of both Coulomb collision and stochastic accel- 
eration, which is valid at sub-relativistic and relativistic 
energies (see e.g. Dogiel (2000) Liang et al. (2002) and 
Dogiel et al. (2007) [) 



df e _ ]_d_f . d/e 

dt p 2 dp \ dp 



B( P )fe 







Here 



B(p) = p 2 



and 



A(p) = B(p) 



( 7 2 _ 1)1/2 



1/2 



P 2 D p (p) 



(24) 



(25) 



(26) 



where D(p) = ap q is the diffusion coefficient due to the 



stochastic acceleration (see e.g. Dogiel et al. (2007) The 
rate of ionization losses is 




7 



In 



kT e ( 7 2 
/ E(p)m e c 2 ('j 2 — 1] 



h 2 uj 2 ^ 2 



- 1)V2 

0.43 



(27) 



where 7 = l + E(p)/m e c 2 7 and w p is the plasma frequency. 
Finally, the electron distribution function takes the form 
(Dogiel et al. 2007) 



f e = C ■ exp ( - 

where 

exp(- 



* B(v) 
— — — - 

A{v) 



Gip) , 



C 



and 



l/0eh 



K 2 (i/e e ) 



(28) 



(29) 



G(p) = 1 - 



Io^xp(J Q w dtB(t)/A(t)) 
jr^™P(fodtB(t)/A(t)) 



The integral 
f P B(v) 
Jo 



A(v) 



-dv 



vdv 




Fig. 5. The spectrum of stochastically accelerated elec- 
trons with acceleration parameter a = 0.0008 (dashed 
curve) is compared with the Maxwcllian spectrum (solid 
curve). A background plasma temperature of kT e = 8.2 
kcV is assumed here. 



The resulting spectrum of the background and accelerated 
particles is shown in Fig. 5. The following limiting cases 
hold: 

i) if a* = 0, the spectrum is a relativistic Maxwcllian 
spectrum 



Yf3 2 exp(--//9 e ) 

e e K 2 (i/e e ) 



(34) 



where K 2 is the modified Bcsscl function of second kind; 
ii) if p <C 1/ s/62, the electron spectrum takes the form of 
the non-rclativistic spectrum derived by Gurevich (1960) 



fe( P ) 



exp 



dw 



-exp 



1/w + a 
dw 



..6 e wi+ 2 



1/w 



1,9+2 



(35) 



(30) 

The position of the cut-off in the electron spectrum is, in 
principle, a free parameter to be fixed at any energy above 
the range of available X-ray observations, e.g. at E > 80 
keV (i.e. the upper value of the energy range of the Beppo- 
(31^AX PDS instrument). In our analysis, we consider specif- 



i+e e v' 



depends on the acceleration parameter 
/£(p)™ e c 2 (7 2 - 1) 



In 



h 2 



wli 2 



0.43 



■) 3/2 j 



(32) 



ically a value of the cut-off energy E cuto f / 



: (7-l) 



The velocity distribution of the accelerated electrons is 
given by 



Ii 2 



fJ 2 (\~p 2 Y/ 2 \^/9 e (l~(3 2 ) 



(33) 



= 212 keV, corresponding to a value of p cu toff = 
where 7 = yl + e p 2 that is able to reproduce both the 
soft and hard X-ray spectrum of Coma (see Dogiel et al. 
2007). In such a model for the origin of the hard X-ray 
emission from the Coma cluster, we derived a value of the 
acceleration parameter a — 0.0008, given the value of the 
electron temperature of Coma kT e = 8.2 keV (Dogiel et 
al. 2007). These parameters fix the shape of the spectrum 
of accelerated electrons and, therefore, the spectral shape 
of the relative SZE. 
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Fig. 6. The dependence of the crossover frequency Xq on Fig. 7. The dependence of the slope S of the SZE in Coma 



the acceleration parameter a given in units of 10 



on the acceleration parameter a given in units of 10 4 . 



4.1. The dependence of Xq and of S on the 
acceleration parameter 

We calculate the dependence on the value of a of both 
Xq and S for the total SZE of the considered cluster 
for the specific model previously discussed. For the orig- 
inal Maxwcllian spectrum of Coma with kT e = 8.2 kcV, 
we obtain values Xq = 3.9 and S = 0.058 in the range 
x = 3.5-4.5. 

We use the relativistic spectrum of thermal electron pop- 
ulation (see Sect. 2) with values of the acceleration pa- 
rameter a in the range a = — 0.0015, bracketing the 
value a = 0.0008 obtained from the fit to the hard X-ray 
spectrum of Coma (Dogiel et al. 2007). 

The approximated expression of the crossover fre- 
quency Xo(a) up to order 0(a 2 ) is given by 

X {a; kT e = 8.2keV) = 3.9 + 15.73a + 15402a 2 . (36) 

For the value a = 0.0008, one obtains X = 3.922. This 
value of Xq would correspond to a thermal electron spec- 
trum (without acceleration) at an effective temperature of 
kT e = 10.5 kcV. 

We note that in this case the shifts induced on the value 
of Xq due to both the electron acceleration process and 
the presence of a peculiar velocity of the cluster are inter- 
twined. The expression for the crossover frequency, taking 
into account the kinematic SZE, is given by 

V- / \ Y / \ _L Vr h i X o( a )] fo 7 \ 

X 0i totai{a) = X (a) H — . (37) 

The peculiar velocity of the Coma cluster was estimated 
to be V r = -29±299fcm/s (Collcss et al. 2001), therefore 
Xo,totai is in the range 3.84 — 4.00 for a value a = 8- 10~ 4 . 
Since the values X^(a = 0) and X^(a = 8 • 10~ 4 ) lie in 
this range, the ability to study the nature of a possible 
population of supra-thermal electrons directly from the 
displacement of the crossover frequency is quite limited. 

The value of the SZE slope (that, we recall, does not 
depend on the peculiar velocity in the frequency range 



x = 3.5—4.5), is therefore an important quantity for deriv- 
ing useful information about the nature of supra-thermal 
electrons. 

The approximated expression of the quantity S(a) up to 
0(a 2 ) is given by 

S(a; kT e = 8.2keV) = 0.0583 + 7.35 • a + 5064 • a 2 (38) 

and its value is 0.0674 for an electron spectrum with 
a = 0.0008. The value of S increases by 16% with 
respect to the case of the Maxwellian spectrum (i.e. with 
a = 0). This value of S would correspond to a genuine 
thermal electron spectrum (i.e. without stochastic accel- 
eration) with an effective temperature of kT = 9.3 keV. 
The dependence of the slope S on the acceleration param- 
eter is shown in Fig. [71 

4.2. Contribution of the kinematic SZE to the slope of 
the total SZE 

The total SZE intensity in the case of a cluster with an 
in-situ accelerated electron population moving along the 
line of sight with a peculiar velocity V r is due to the com- 
bination of SZth and SZkin 

i(x) = ith(x) - —h(x) . (39) 
c 

The series expansion of the term SZ t h around the adimen- 
sional frequency x = 4 is given by 

ith = km + S-{x-A) + 0{x- 4) 2 , (40) 

where S is the value of the slope of the SZE without the 
contribution of the kinematic SZ effect. 
The series expansion of h(x) around x = 4 is 

h{x) = 4.865-0.181(x-4)-0.512(x-4) 2 +O(2;-4) 3 ,(41) 

which can be written as 

h{x) = h + A(x) ■ (x - 4) , (42) 
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where A(x) = -0.181 - 0.512(a; - 4) + 0{x - 4) 2 . 
The total SZE intensity is then given by the expression 



i{x) = itho 



— ho 

c 



For a narrow frequency range, the spectral slope of the 
SZE can be derived by means of SZE measurements at 
two nearby frequencies, given by 



S 



bs 



i(xi) - i(x 2 ) 



(44) 



Xi - x 2 

For the sake of clarity, we fix the first frequency X\ = 4 
and the second frequency lies in the range x 2 = [3.5,4.5]. 
In this case, we find that 



Sobs = Sth - — A{x 2 ) 
c 




-0.0015 



Fig. 8. The contribution of the kinematic SZE to the slope 
of the total SZE (see text for details). 

The contribution of the kinematic SZE to the slope of the 
total SZE is given by 



Skin S t 



obs 



S, 



Hi 



10 3 km/s 



AS{x 2 ) 



(46) 



where the quantity AS(x 2 ) = 3.3 • 10 3 A(x 2 ) is plotted 
in Fig. E 

It is clear that the contribution of the kinematic SZE to 
the slope of the total SZE is far smaller than the slope 
increase due to the stochastic acceleration process (see 
Fig. [7]). It follows that the determination of the slope of 
the SZE is a powerful tool for studying the presence and 
the nature of a possible population of accelerated electrons 
at supra-thermal energies in the atmospheres of galaxy 
clusters. 

5. Discussion and conclusions 

We demonstrated in this paper that the value of the 
crossover frequency Xq of the SZE depends crucially on 
the cluster peculiar velocity V r (as well as on the CMB 



primary anisotropies), while the value of the slope of the 
SZE does not depend on the kinematic SZE spectrum (and 
on the CMB anisotropy spectrum), in the optimal fre- 
quency range around the crossover frequency of the ther- 
mal SZE, i.e. in the frequency range x = 3.5 — 4.5. This 
is because both the spectrum of the kinetic SZE and that 
of the CMB primary anisotropies are almost flat in this 
frequency range. Therefore, while the amplitude of the 
SZfej„ produces a systematic bias in the position of the 
crossover frequency Xq, it does not affect significantly the 
slope S of the SZE. 

This fact is more evident when one computes the CMB 
temperature change produced by the SZE, i.e. the quan- 
tity 



(45) -=- 



AT _ (e x - l) 2 AI 



To x A e x 



(47) 



in fact, the CMB temperature change induced by the ki- 
netic SZE has a flat spectrum in frequency, while any 
other CMB temperature change induced by other SZE 
have spectra with specific values of slopes that differ from 
zero. 

Since the cluster peculiar velocity is unknown in many 
cases, the precise position of the crossover frequency Xq 
of the thermal effect cannot be well-defined from the SZE 
observations because the measurements of peculiar veloc- 
ity are highly uncertain. 

To overcome this systematic bias, we have proposed 
to use the spectral slope of the SZE to obtain unbiased 
information about specific properties of various electron 
distributions in galaxy clusters: the gas temperature for a 
thermal electron population, the power-law spectrum of a 
non-thermal electron population, the combined spectrum 
of thermal and non-thermal electron populations, the ac- 
celeration parameter of a stochastically accelerated elec- 
tron population. We have evaluated analytically specific 
values of the slope for various SZE signals in the frequency 
range x = 3.5 — 4.5. We found that in this frequency range 
(around the crossover frequency of the thermal SZE) the 
spectral slope of the SZE does not depend on the value of 
the cluster peculiar velocity because the spectrum of the 
kinetic SZE is extremely flat in this frequency region. We 
have shown this for different electron populations residing 
in clusters. 

A doubtless advantage of the method that we have 
presented here is that it provides an opportunity to search 
more reliably for nonthermal components of electron spec- 
tra in clusters and determine their characteristics. These 
allow us to derive independently parameters of accelera- 
tion processes in clusters which are usually unknown and 
estimated in the framework of different models. This fi- 
nal circumstance is particularly important for the problem 
of the origin of nonthermal emission from clusters in the 
Extreme UV, in the hard X-ray and at radio wavelengthes. 

Spectroscopic measurements of the SZE spectra will 
become available in the near future by means of 
dedicated space experiments (such as e.g. SAGACE; 
Spectroscopic Active Galaxy And Cluster Explorer, see 
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http:/ /oberon. romal.infn.it/sagacc/) based on microwave 



spatially resolved spectroscopic techniques. The advent of 
these spectroscopic capabilities will allow us to study in 
an unbiased way several physical details of the leptonic 
structure of the atmospheres of galaxy clusters and other 
cosmic structures. 
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Appendix A: The shape of the SZE around the 
crossover frequency 

In this Appendix, we calculate the values of the crossover 
frequency Xq and the slope S of the SZE in the frequency 
range x = 3.5 — 4.5. We evaluate these quantities for vari- 
ous cases of thermal and non-thermal electron populations 
with the inclusion of the kinematic SZE. 



A.l. Thermal electron spectrum 

The total SZE intensity, in the case of a cluster with a sin- 
gle thermal electronic population moving along the line of 
sight with a peculiar velocity V r is due to the combination 
of SZ t h and SZfci„ and is given by 



i(x) 



m e c 2 



V r m e c 2 . 



In the non-relativistic case, the series expansion of g(x) 
and of h(x) around x = 4 (i.e. the central value of the 
frequency range x ~ 3.5 — 4.5) is given by 

g{x) = 0.726 + 4.28(2; - 4) - 0.038(2; - 4) 2 - 0.489(2; - 4) 3 
+0(x - 4) 4 (A.2) 



and 

h(x) = 4.865 - 0.181( 
+0{x-4) 4 



(x - 4) - 0.512(;r - 4) 2 + 0.088(:c - 4) 3 

(A.S) 



Thus, the series expansion of i(x) around x = 4 is given 

by 



i(x) 



kT„ 



0.726-4.865 



V r m e c' L 
c kT e 



4.28 

kT e 
m e <y 

4.28 



3 l^I^] (a ._4) 



c kT ( 
0.726- 0.827 



V r lOkeV 



0.0308 



10 3 fcm/s kT c 
V r lOkeV 



Wkm/s kT e 



(x-4) 



,(A.4) 



where we neglect higher order terms in the frequency range 
x ps 3.5 - 4.5 (i.e. v ps 200 - 255 GHz). In the first brack- 
ets of Eq. (| A.4|) , the term depending on V r and T e cannot 
be neglected, but, in the second brackets of the same Eq. 
(|A.4j) we can neglect the term depending on the combi- 
nation V r m e c 2 /ckT e , since it is always far smaller than 
the value 4.28 for realistic values of V r and kT e found for 
clusters. 

The slope, S= [i{x)-i(X )]/(x-X ) =i(x)/(x-X ) 
of the non-relativistic SZE spectrum is 



kT e g{X$ + Sx) 
m e c 2 Sx 



nf-L 'f kTi 

4.25 -7T = 0.08 



m e c^ 



lOkeV 



(A.5) 



Hence, the slope of the SZE intensity i(x) in this frequency 
range does not depend on the cluster peculiar velocity. 

The value of Xq of the total (thermal plus kinematic) 
SZE is defined by the condition i(Xq) = 0, or 



kT e 



0.726-4.865- 



. V r m P c 



+ 4.28(X -4) 



From the previous condition, one obtains 
0.726\ 4.865 V r m e c 2 



X 



o.total.nr 



4.28 J ' 4.28 c kT e 
which is written, in more practical units, as 

V r lOkeV 



= 0.(A.6) 



(A.7) 



X 



o.total.nr 



3.83 + 0.193 



10 3 kms 



kT„ 



(A.8) 



In the rclativistic treatment, the CMB spectral distor- 
tion due to the thermal SZE is given by the expression 



(A.l) Alt 



n 1 



2(kT Q f 
h 2 c 2 



dsPi(s) 



cxp(2;e 



l(A.9) 
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where 



Pi(«)= / d(3f e (P)P(s,/3) 



(A.10) 



with f e (P) being the velocity spectrum of the electronic 
population and P(s, (3)ds being the probability that a sin- 
gle scattering of a CMB photon off an electron with speed 
[3c causes a logarithmic frequency shift s = ln(y' jv). Here, 
we use the first order approximation in r of the pho- 
ton redistribution function P(s) for the sake of simplic- 
ity. Analogous conclusions hold, nonetheless, for the more 



general derivation (see e.g. Colafrancesco et al. (2003) 



The series expansion of i{x) around the crossover fre- 
quency Xq is written as 

XQexp(Xoexp(— s))exp(— As) 



i{x) = / dsPi(s) 



(exp(Xoexp(— s)) — l) 2 
X$exp(X ) 3 / X^exp(-3s) 



(exp(Xo) — l) 2 Xq \ exp(Xoexp(—s)) — 1 

H(x-Jr ). (A.11) 

exp(X )-l)\ { 0> V ' 

Using th e expression X (T) = 3.830(1 + 1.1620 e - 
0.81446> 2 ) [Dolgov et al. (2001| one can calculate the slope 
of the thermal SZE for the relativistic treatment as a func- 
tion of kT e . 

A good approximation of the slope for the relativistic 
treatment, and in the temperature range kT e = 1 — 20 
keV, is given by 



Jrel 



4.25- 



0.08- 



kT e 
m e c 2 
kT e 
WkeV 



30 



kT e _ 



1 - 0.138 



kT e 
\0keV 



(A.12) 



The value of crossover frequency of the total SZE in the 
relativistic treatment X ^totai,rei is given by 



,total,rel ~ 

X (T e ) + 

« X (T e ) + 



V r h(X (T e )) 
C Srel 

0.193 



1-0.138- 



fcT e 
WkeV 



V r lOfceU 



Wkms- 1 kT e 

A. 2. Power-law, non-thermal electron spectrum 

For a power-law, electron spectrum given by 



(A.13) 



fe( P )= l-r\-. -P~ a Pl<P<P2, 



hip) 



o 



elsewhere 



(A.14) 



with p = /?7, the electron spectrum in terms of f3 is given 
by 



3/2 



(A.15) 



If p\ is ultra-relativistic (i.e. p\ 3> 1), the limiting form of 
i(x) is given by 



^ultra — rel \Xj 



(A.16) 



e x - 1 

One can derive analytically the approximate dependence 
of Xq on 7 for a monoenergetic electron spectrum. The 
photon density distribution for a CMB-blackbody spec- 
trum is 

e x — 1 

An electron with Lorentz factor 7 increases the frequency 
a; of a scattered photon on average by a factor x'/x = 
47 2 /3. For relativistic electrons, the photon is therefore 
scattered to much higher energies. The photon spectrum 
after a scattering is 



n'(x') = n[a;(x')] 



dx 
dx 1 



1 



3x' 

47 2 J exp(3x'/4:j 2 ^ 



3 

T ' i^ 2 



(A.18) 



The value of 3x'/47 2 is <§; 1 for 7 ^> 1. Since one expects 
in this case that Xq is much larger than 1 because of large 
amounts of transferred electron momenta, we can find the 
value of Xq from the condition that exp(x' ) 3> 1. In this 
case, the value of Xq is defined by the condition 

K) 2 



n'(x' ) - n(x' ) 
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I67 



4^0 



exp(x' ) 







(A.19) 



Since Xq 1, the condition X 3> ln(Xo) holds, and one 
obtains 



X a ^2ln[ - 7 2 



(A.20) 



A.3. Combination of thermal and non-thermal electron 
populations 

We evaluate the values of Xq and S for a combination 
of a thermal plus a power-law electron spectra with the 
addition of a peculiar velocity of the cluster. 

One can consider a photon to be removed effectively 
from the CMB spectrum if its energy is increased by about 
one order of magnitude or more, requiring 7 > 3 (e.g. 



Ensslin & Kaiser (2000) Using the usual Kompaneets ap- 



proximation and neglecting those relativistic electrons 
with 7 < 3, the produced distorted spectrum is 



i{x) 



kT P , 



h(x)- 



7~rel fl1 e C X 



kT e t kT e e x — 1 

The series expansion of i{x) around x — 4 is then 



(A.21) 



i{x) 



kT e . 



0.726-4.865- 



V r m e c 



4.28 + 0.181 



V r m e c 
c kT e 



kTe 



0.321 



1.194- 



Trel m e C 

~~kT, 



rel TTLeC 
T kT e 

>-4) 



(A.22) 
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From the condition i(x) = 0, the value of the crossover 
frequency X of the total SZE is 

1 0.726 - 4.865-^=^! - 1.194^-^ 

X ° = 4 ~4~28 X 1 i 0-181 V r m e c 2 , 0.321 r rcl 7n e c~ (A. 23) 
1 ' 4.28 c fcT„ ' 4.28 r fcT e 

From Eqs. (|A.22[) and (|A.23|) . we then derive the values of 
0.726 4.865 V r m e c 2 1.194 T rel m e c 2 



X = 4 



4.28 4.28 c kT e 4.28 r /cT e 
7 r 10fceV s 



3.83 + 0.193 



10 3 /cm/s kT e 



O.U8[ Trel/TWk f) (A.24) 



0.01 fcT e 
and of the slope of the total SZE 

5 = 4.28^1 + 7.10-3 * WkeV 



m e c 2 \ 10 3 km/s kT e 

Trel It lOfeel 
0.01 kfl 



038 Z^/r 10^F\ (A 25) 



